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Abstract 

We study a general problem of the translational/rotational/vibrational/electronic dynamics of a diatomic 
molecule exposed to an interaction with an arbitrary external electromagnetic field. The theory developed 
in this paper is relevant to a variety of specific applications. Such as, alignment or orientation of molecules 
by lasers, trapping of ultracold molecules in optical traps, molecular optics and interferometry, rovibrational 

>: 

spectroscopy of molecules in the presence of intense laser light, or generation of high order harmonics 
Q\ ■ 

^— < from molecules. Starting from the first quantum mechanical principles, we derive an appropriate molecular 

7-h ; 

Hamiltonian suitable for description of the center of mass, rotational, vibrational and electronic molecular 

motions driven by the field within the electric dipole approximation. Consequently, the concept of the 

Q_iBorn-Oppenheimer separation between the electronic and the nuclear degrees of freedom in the presence of 
i 

£3 an electromagnetic field is introduced. Special cases of the dc/ac field limits are then discussed separately. 

a \ 

^^Finally, we consider a perturbative regime of a weak dc/ac field, and obtain simple analytic formulas for the 
j> associated Born-Oppenheimer translational/rotational/vibrational molecular Hamiltonian. 
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I. INTRODUCTION 



During the last decade the manipulation of molecules by lasers has been extensively studied 
both theoretically and experimentally. The alignment and orientation of molecules by lasers have a 
large variety of applications in different fields of chemistry, physics, and potentially also in biology 
and material research [jj. Examples of recently demonstrated applications range from laser- assisted 
isotope separation 0] and catalysis Q|, from pulse compression ^ and nanoscale design 0, 0] to 
tomographic imaging of molecules [t| and quantum information processing jsj]. 

The Hamiltonian for molecules in an external electromagnetic field is of interest since it underlies 
a variety of phenomena associated with the electromagnetic field control of external and internal 
molecular motions, including trapping [j]], molecular optics ID, JjlQ. Stark shift manipulation 
of the potential energy surfaces 3], and control of the high order harmonic generation 

Surprisingly, two qualitatively different forms of the rovibrational Hamiltonian for molecules in 
weak laser fields appear in the theoretical literature dealing with laser alignment. One form of the 
Hamiltonian has been derived in Ref . 3] and the other one in Ref . . Both approaches, although 
being contradictory, have been used extensively in theoretical studies, giving thus rise to serious 
confusions and controversies. Most recently we have resolved this "puzzle" ^3] by applying the 
adiabatic theorem for open systems using an extension of the {t,t') method, termed the (t,t',t") 
approach 

The purpose of this work is to provide a detailed derivation of the Hamiltonian for diatomic 
molecules in laser fields, regardless if the involved field intensities are weak or strong. Our motivation 
is to analyze the most general case of a " diatomic molecule - electromagnetic field" interaction, such 
that the obtained results should be relevant not only for description of molecular alignment, but also 
for trapping of cold molecules in optical lattices, for rovibrational spectroscopy of molecules in the 
presence of intense laser light, or for generation of high order harmonics from molecules. 

The paper is organized as follows. In Section II, we present a rigorous self contained derivation 
of an appropriate molecular Hamiltonian suitable for description of the center of mass, rotational, 
vibrational and electronic molecular motions driven by the field within the electric dipole approxima- 
tion. Consequently, in Section III we introduce the framework of the Born-Oppenheimer separation 
between the electronic and the nuclear degrees of freedom in the presence of an electromagnetic field. 
Concept of a time dependent electronic potential energy surface is then discussed, with particular em- 
phasis on the special cases of the dc/ac field limits. In Section IV, we consider a perturbative regime 
of a weak dc/ac field, and establish an interconnection between the corresponding Born-Oppenheimer 
electronic potential energy surfaces and the conventionally used static /dynamic molecular polariz- 
abilities. Concluding remarks are given in Section V. 



II. DIATOMIC MOLECULE IN AN ELECTROMAGNETIC FIELD: 
THE HAMILTONIAN 



A. The Hamiltonian in momentum gauge and in laboratory frame coordinates 



Let us study a diatomic molecule AB exposed to an interaction with laser light. Some external 
electrostatic field can also be present. We prefer here to describe the considered electromagnetic field 
classically, in terms of the scalar potential <j)(r, t) and the vector potential A(r, t) , using Coulomb 
gauge and Gaussian units for the electromagnetic quantities The corresponding molecular 

Hamiltonian (expressed with respect to the laboratory space fixed coordinate frame) possesses an 
explicit form 
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Here, symbol e stands for a charge of an electron, c denotes the velocity of light, Z A and Z B are 
the atomic numbers of the two nuclei A and B, while terms m A ,m B and m e represent respectively 
the masses of nuclei A, B or the mass of an electron. An auxiliary index j = 1,2, ... ,N has been 
adopted for labelling the electronic variables. Other notations should be self explanatory. 



B. The momentum gauge Hamiltonian in the center of mass and relative coordinates 

As the first step of our derivation, we switch from the laboratory frame coordinates into the center 
of mass and relative coordinates. To accomplish this task, we introduce the center of mass position 
vector 

3 m A f A + m B f B + J2f=i m e r j 

■tic = TT i l^J 

M y ' 

where the total mass 

M = m A + m B + N m e . (3) 
In addition, we define the relative coordinates 

Rab = r A - f B ; (4) 



and 



Relations inverse to the formulas (J2J), (HJ) and (J3J) are easily found to be 
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and 

r} = R c + % . (8) 
An additional auxiliary symbol has been adopted here, 

m AB = m A + m B . (9) 

For the sake of completeness, we also mention in the present context that the volume element remains 
unchanged after the above described coordinate transformation, i.e., 

N N 

d 3 r A d 3 r B Y[d 3 rj = d 3 R c d 3 R AB \{d 3 q j . (10) 
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Proceeding further, we introduce the momenta associated with the new coordinates. Namely, we 

define the operators 

p c = -ihv nc , P AB = -iftVa AB , pj = . (n) 

These new momenta are interconnected with the original laboratory frame momenta through the 
transformation formulas 
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and 
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P J = W Pc + ft ~ W ^ ' (14) 
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We continue by substituting Eqs. ©, ©, © and Eqs. (fI2jl . (IT3|) . (fill into Eq. (0). In order to 
simplify the obtained result, we employ the dipole approximation 

A(R c + lt) « A(R c ,t) , </>(Rc + U) « 0(/M + g-V Mc <f)(R c ,t) ; (15) 

which is justified as long as the spatial variation of the electromagnetic field remains negligible at 
the length scales |£| comparable to molecular dimensions. By using also the transversal property of 
the Coulomb gauge vector potential 

V-A(f,t) = ; (16) 
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we write down an explicit expression for the Hamiltonian, Eq. (0), in the center of mass and relative 
coordinates. It holds 
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In Eq. (117)1 an auxiliary shorthand symbol 



m A m B 



m A + m B 

stands for the reduced mass of the AB molecule. Additional simplifications are in order: i) The term 
(l/(2M))J2jj' Pj " Pj' in Eq. (JT7J) can be neglected since the factor (1/M) is small in magnitude. 
Similar argument applies also in the case of terms (m e /m AB ) Y^jLi Qj ■ ii) F° r neutral molecules the 
[Z A + Z B — iV]-dependent contributions to Eq. (11711 v anish. Hi) The A 2 (R C , t) factor can be eliminated 
from Eq. (|17|1 by a trivial phase transformation |20j], provided that we neglect additional corrections 
of the form M~ x [ /* A 2 (R C , t') dtf ] ■ P c and M- x [ A^f A 2 (R c ,t') dt') which arise due to non- 
commutativity between A(R c ,t) and P c . This step is justified as long as the spatial derivatives of 
the vector potential remain sufficiently small such that the translational motion of the molecule is 
not affected by the mentioned correction terms. 

Having incorporated all the above simplifications, we rewrite the Hamiltonian ()17j) into a relatively 
simple functional form 
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The center of mass motion becomes here nonseparable from the internal molecular motions solely 
due to presence of the field terms A(R C , t) and <f)(R c ,t) in Eq. (fT^j). 



C. The length gauge Hamiltonian in the center of mass and relative coordinates 



As the second step of our derivation, we convert the Hamiltonian (|19|) into the length gauge 
which lends itself better for practical applications discussed later in Section III. The length gauge 
Hamiltonian H(t) is obtained by an unitary transformation 



H(t) = U+(t) H(t) U(t) - ihU\t)^U(t) 
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with an unitary operator 
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Straightforward algebraic manipulations reveal that 
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Here, the quantity 
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can be interpreted as the dipole moment operator of AB molecule, and symbols 
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c dt 

stand for the transverse and the longitudinal electric fields assigned to the potentials A(R C , t) and 
<f)(R c ,t) , respectively [3]. For the sake of completeness, we note by passing that in the formula 
we have actually neglected additional corrections arising due to non-commutativity between 



the operators A(R c ,t) and . Justification of this step is the same as in item in) of the previous 
subsection II. B. 

Before proceeding further in our derivation, let us mention a few interesting observations regarding 
the quantity (|23|). For homonuclear molecules (m_4 = mg and Z A = Z B ) the first term of equation 
(|23j) vanishes and thus only the electronic contribution [— e Y^f=i Qj] is relevant. On the other hand, for 
cases where = Z B = Z but ^ mg due to the use of different isotopes (such as HD for example), 
the formula (J23J) contains a factor Ze (mg — mj^/m^AB R^b which is acting as a "permanent-like" 
dipole moment and influences the photo-induced molecular dynamics. The mentioned "permanent- 
like" dipole moment contribution arises in the case of isotopically substituted homonuclear molecules 
solely due to the fact that the nuclear center of mass is not located in the geometrical center of the A— 
B bond (which constitutes a molecular symmetry center from the point of view of electronic structure 
calculations). One might expect that the above discussed dipole moment component e [(2^mg — 
l Z'B m A) l m AB\ Rab becomes even more important in the case of heteronuclear diatomics (m_4 7^ m# 
and Z A ^ Z B ). 



D. The length gauge Hamiltonian in the spherical polar coordinates 



As the third step of our derivation, we replace the three cartesian coordinates Rab — 
(Xab, Yab, Zab) by their spherical polar counterparts (R, $,<p). We employ the usual transfor- 
mation procedure which is well known e.g. from standard textbook treatments of the hydrogen atom 



2l| . The corresponding transformation formula reads as 
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The associated volume element is of course d 3 R AB = R 2 dR sin & dd df . What remains to be done 
is to rewrite the Hamiltonian (|22|) into the new coordinates. An appropriate procedure for resolving 



this task is well established, see e.g. Chapter IX of Ref. |2j. Therefore, we display here explicitly 
just the final result, 
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Here, the radial momentum operator is defined as 
the squared angular momentum operator is given by 

h 2 



i \qj + Nm^)M(^)fl5s 



f 2 



sin 2 19 

and an additional auxiliary symbol 
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To avoid confusion, let us note explicitly that the dipole moment operator (j2SJ) is now expressed in 
the form 
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E. The length gauge Hamiltonian in the body fixed electronic coordinates 

In this step of our derivation, we transform the position vectors of all the electrons into the 
body fixed frame. The origin O of the body fixed coordinate system is set to be the molecular 
center of mass. Note that this choice of the origin is a bit different from the choice adopted within 
the usual spectroscopic literature, where the nuclear center of mass is considered instead (see for 
example Ref. We prefer to use here an alternative less conventional assignment of the origin 

O since it makes our formulation more transparent and enables us to avoid introducing additional 
appr oximat ions . 

The body fixed © o z axis is, by definition, parallel (although not always coincidental) with the 
direction of R AB . The body fixed o x and o y axes are constrained by the requirement © o x x 
© o y = © o z . Choice of o x and o y is, however, not unique: An arbitrary rotation around o z leads 
to an equivalent pair of body fixed axes (o' x ,o' y ) which are equally suitable as (o x ,o y ) . Hence, an 
unambiguous definition of o x and o y must be fixed by convention. In order to achieve maximum 
simplicity, we prefer to employ such a particular convention that 



q 3 - = M{ti,(p) Tj 



(33) 



where fj are the body fixed coordinates of vector <jj . Since the matrix Ai($,ip) is orthogonal, 
the volume element remains unaffected, d 3 qj = d 3 Tj . Having introduced the body fixed electronic 
coordinates, we continue further and define the associated momenta, 



?j = -ihVtj . (34) 

These new momenta are interconnected with their space fixed counterparts through the transforma- 
tion formulas 

Pj = M(-&,<p)vj • (35) 

It is straightforward to rewrite the Hamiltonian (f2"K|) into the body fixed coordinates. Taking advan- 
tage of the orthogonality property (J27)) . we arrive towards the desired result 

2 f*2 N -2 ry V O 



"P 2 I) T I N ^2 7 7 

U " 2M 2fi AB 2fi AB & g2m e R r s r> 



e 



J=i \/ x J 2 + y i + t z J ~ ( m e/^B)^] 2 £i yjxj + y] + [zj + {rn A /m AB )R} 2 
- \M{$,ip)D%{R,$,ip^ N )]-{E\R c ,t) + E\R c ,t)} . (36) 

Here, the quantity 
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represents the body fixed counterpart of the dipole moment vector ()23|) . Note that in formula (f3l)j) the 
translational and rovibrational kinetic energy operators are completely decoupled from the kinetic 
energy operators of the electrons. This holds true in particular also for the electronic and the nuclear 
angular momenta. (The electronic angular momenta are not displayed here explicitly and appear 
only after switching into the spherical or cylindrical electronic coordinates.) 



F. Final form of the Hamiltonian for diatomic molecules in an electromagnetic field 

Summarizing all the elaborations of Section II, we may conclude that the quantum dynamics of 
the considered molecule AB interacting with an external electromagnetic field A(f, t) and 4>(f,t) is 
described by the time dependent Schrodinger equation 

ih ^ E(R C: R, 0, <p, r N : t) = H(t) E(R C , R, 0, <p, r N : t) ; (38) 

where E(R C , R, 0, if, r N , t) is the associated translational/rotational/vibrational/electronic wavefunc- 
tion, and the appropriate Hamiltonian H(i) is given by expression ()36|) . To avoid confusion, let us 
note in passing that the electron spin variables are suppressed in the notation of the present paper, 

n 



since they never enter explicitly into our considerations. Nevertheless, the presence of an electronic 
spin is of course fully respected within our treatment, as well as the antisymmetry of the electronic 
wavef unctions. 

Before proceeding further, it is convenient to introduce an additional simplification, based upon 
the factorization 

5(4,^,0, = R-^iRcR^^^t) ■ (39) 

The purpose of this factorization is to eliminate redundant difficulties arising due to a compli- 
cated functional form of the radial momentum (|29j). We refer again to standard textbooks 2jJ 
for a more detailed discussion of this issue. One can easily show that the redefined wavefunction 
^(R c , R, <f, r N , t) satisfies the time dependent Schrodinger equation 
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with the Hamiltonian 
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Here, the field free electronic Hamiltonian 
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the " molecule - field" interaction term 

W el (R c ,R,&,<p,t) = - [M{d,<p)D%%{R,#,<p,T N )]- E(R c ,t) ; (43) 
and an overall electric field 

E(R c ,t) = E\R c ,t) + E\R c ,t) . (44) 



III. DIATOMIC MOLECULE IN AN ELECTROMAGNETIC FIELD: 

THE TIME-DEPENDENT BORN-OPPENHEIMER ELECTRONIC POTENTIAL 

ENERGY SURFACES 

A. The time-dependent electronic wavefunctions 

For isolated molecules, the well known concept of the Born-Oppenheimer/adiabatic separation 
between the electronic and the nuclear degrees of freedom proved to be extremely useful, as it gives 



a lot of physical insight by distinguishing between the electronic and the rovibrational molecular 
states. The purpose of this Section is to extend the formulation of the Born-Oppenheimer approach 
on cases when a molecule is exposed to an interaction with an electromagnetic field. For simplicity, 
we shall consider here just a diatomic molecule AB for which an appropriate Hamiltonian has been 
discussed at length in the previous Section II. 

As the first step of our analysis, we formally construct the time dependent Born-Oppenheimer 
electronic basis set, defined by particular solutions of the electronic time dependent Schrodinger 
equation 

ih^- t $ n (r N ,t;R c ,R,#,tp) = [H e i(i?) + W e i(.R c , R, i?, tp, t) j $ n (T N , t; R c , R, i?, ip) . (45) 
The corresponding initial condition is conveniently chosen to be 

<f> n (r N ,t ;Rc,R,'&,<p) = K(? N ;R) ; (46) 

where t is an as yet arbitrary time instant, and function $° (r N ; R) represents a solution of a field 
free electronic eigenproblem 

K el (R)<S>° n (r N ;R) = S° n (R) ^ n (r N ; R) . (47) 

In order to simplify the underlying notation, we assume here that the electronic Hamiltonian H i(i2) 
possesses a discrete spectrum labelled by a single collective index n . Note, however, that the con- 
tinuum part of the spectrum of H e i(i?) can be implicitly included in this way as well, by taking 
advantage of the box quantization procedure. For the sake of clarity, let us also recall once again 
that the electron spin variables have been suppressed in above formulas, although they are implicitly 
accounted for. 

The wavefunctions defined by Eq. (|4"5j) depend only parametrically on the nuclear coordinates 
(R c , R, ip) . For each fixed nuclear configuration (R c , R, <p) , the associated collection of electronic 
states {|$ n (t))} forms a complete orthonormal basis set covering an entire Hilbert space of the 
electronic variables. The orthonormality and closure properties are granted here for every time 
instant t, since the time evolution according to the Schrodinger equation, Eq. (|45jl . is unitary, and 
since the field free electronic eigenf unctions (JTTj) entering into the initial condition ()46|) constitute 
themselves a complete orthonormal electronic basis set. 

We note by passing that although the above initial condition ()46|) can be in principle used in the 
most general context, it carries an especially physically illuminating interpretation in the case when 

A(R c ,t) = , <p(Rc,t) = , W el (R c ,R,$,<p,t) = (t < t ) ; (48) 

valid for all the possible molecular positions R c under study. If so, each stationary quantum state of 
the considered molecule can be (of course for t < to and in the absence of avoided crossings) char- 
acterized within the framework of the conventional Born-Oppenheimer approximation, as a product 



of a specific electronic state |$°) exp[— {i/fi)£ ®(R)(t — t )] and an appropriate nuclear component. 
When the field is switched on at t > t , it is natural to expect that the relevant electronic state n 
remains relatively well defined, provided that its time evolution is understood in the sense of equation 
(|45jl . This idea stands behind our formulation of the generalized time dependent Born-Oppenheimer 
separation scheme which is elaborated in the coming subsection. 

B. Generalized Born-Oppenheimer theory for time dependent Hamiltonians 

Following the spirit of the Born-Oppenheimer approach, we expand the total molecular wavefunc- 
tion as 

y(R c ,R,$,<p,r N ,t) = J2Xn(Rc,R,#,<P,t)$ n (? N ,t-,R c ,R,&,<p) . (49) 

n 

Substitution of an ansatz (|4"9|l into the time dependent Schrodinger equation (J4"U|) leads towards a 
set of coupled differential equations for the as yet unknown expansion coefficients \ n . That is, 
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where the operator 

h nn ,(t) = h£>,(f) + h« (t) + h2,(t) (51) 
is defined in terms of its components through the formulas 



2M R c 2fi AB dR 2 2fi AB R 2 



(52) 



hS(t) = - ^ ($ n |V^|^> ? , ■ - £- (<S> n \(d/dR)\<S> n/ ) ?N A (53) 
" (*„|(0/0*)l*»'>** ^ ($ n |(9/^)|^) ? - |- ; 

hS(t) = ($„|hS(t)|$ n ,) fJ v = -J^^A^)** - ^-^ n \{d 2 /dR 2 )\^ n ,), N 

fi 2 Ti 2 cn^i9 

($ n \(d 2 /d$ 2 )\$ n ,)? N - -f—- (<S> n \(d/d#)\<S> n ,)t N 



2\i AB R 2 2fi AB R 2 shit? 



2fx AB R 2 sin 2 ?? 

Here, for the sake of notational compactness, the matrix elements over the electronic wavefunctions 
have been expressed in a condensed fashion, such that e.g. 

($ n \(d/dR)\$ n ,)?N = J d 3 r N ®* n (r N ,t-,R c ,R,#,<p) — $ n ,(r N ,t; R c , R,#,cp) ; (55) 

and similarly for the other quantities of this kind appearing in Eqs. (|53| ) -([54 p . 



Relations (|49|) and ()50|) are still exact, as they constitute just a different equivalent representation 
of the original problem ([40)1 . Now, let us introduce the generalized adiabatic approximation, based 
upon neglecting all the off-diagonal (n ^ n') terms h£2,(£) and h^],(£) of the Hamiltonian. If so, 
an index n (used for labelling single dynamical electronic basis vectors) becomes a good quantum 
number, and equations (|49|) and (J5T)|) are simplified into 

vl> BO (^ c , R, 0, <p, r N , t) = Xn°(Rc R, 0, <f, t) d? n (r^ t; R c , R, 0, <p) ; (56) 

and 

ih^Xn°(Rc,R,#,<P,t) = h nn (t)x*°(Rc,R,#,<p,t) ■ (57) 

Before proceeding further, let us briefly touch a question regarding the justification of the above 
approximative step. In analogy to the conventional case of the time independent (field free) Born- 
Oppenheimer approximation, it seems reasonable to suggest the perturbation method as a well 
defined way how to establish a criterion of validity for an adiabatic ansatz l|56j) . If so, one might 
speculate that the adiabatic approximation (|56| ) - (j57|) breaks down in presence of "near-degeneracies" 
whose enhancement or suppression can be strongly influenced by an electromagnetic field. An explicit 
elaboration of the just sketched ideas seems relatively straightforward but lies beyond the scope of 
the present text. 

C. Time dependent electronic potential energy surfaces: The general case 

Equation (JSTj) is naturally interpreted as an effective time dependent Schrodinger equation which 
governs the nuclear motions of the AB molecule in a given dynamical electronic state n. It is conve- 
nient to carry out an additional transformation which converts the associated effective Hamiltonian 
h nn (t) into a more standard functional form understandable as a sum of the kinetic and the potential 
energy operators. To accomplish the mentioned task, one needs to eliminate the Hamiltonian term 
n nn W containing the first order derivatives with respect to the nuclear coordinates. An appropriate 
procedure for this is based upon making a factorization 

X*°(Rc,R,^^t) = e w ^ R ^URc,R,#,<P,t) ; (58) 

where an exponential factor w n (R c , R, 0, ip, t) is chosen to satisfy a system of the first order partial 
differential equations 

V nc w n {R c ,R,#,ip,t) = -<$„|V£j$ n )** , dw n /dR = -{$ n \{d/dR)\$ n )?N (59) 
dw n /d# = -<$„|(<9/dtf)|$ n ) ?J v , dw n /dip = -($ n \(d/dtp)\$ n )?N ; (60) 



valid at every time instant t . 

The question on general solvability of Eqs. (|59| ) -()60 p is nontrivial and it is out of the scope of 
this paper to discuss it. We shall assume for now that the desired solution w n (R c ,R,'&,<p,t) exists 
and is unique up to a constant factor. This is of course the case of a field free problem, where the 
time dependent electronic wavefunctions depend only on R and not on the other nuclear coordinates 
(R c , f) . We refer also to other important special cases discussed in subsequent subsections III.D 
and III.E, where an exact or at least a well defined approximative solution can be shown to exist. 
In this context one should note that, even if the set of Eqs. (|59 |) -(|S0 j) turns out to be generally not 
solvable, an essential concept of the Born-Oppenheimer separation represented by Eqs. (|56)l - (|57) l 
remains unaffected. The only controversy persists here on whether or not one can formulate the 
general Born-Oppenheimer theory in terms of the electronic potential energy surfaces. If not, one 
has to use the above untransformed Born-Oppenheimer Hamiltonian h nn (t) which contains first order 
derivatives with respect to the nuclear coordinates. 

Having fixed the exponential factor w n (R c , R, <p, t) according to Eqs. (|59 Jl -(|S0 ]) . an equation of 
motion for the wavefunction component £, n {Rc, R, f, t) is found to be 



ih?- t URc,R,#,<P,t) = h^(t)URc,R^,^t) 

where an effective translational-rotational-vibrational Hamiltonian 



(61) 



rr 



A, 



h 2 d 2 



+ V n {R c ,R,ti,<p,t) 



(62) 



2M Rc 2/2 AB OR 2 2/j AB R 2 
can be, indeed, understood as the kinetic energy operator of the nuclei plus an electronic potential 



term. An explicit expression for the mentioned potential is 



V n {R c ,R^^,t) 



d ~ .-. 

ih — w n (R c ,R,^,(f,t) + V n (R c ,R,$,tp,t) 
at 



(63) 



with the term 



V n (R c ,R,#,<p,t) = \h^(t) + hU(t)]w n (R c ,R,#,<p,t) + h£>(t) 
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h 2 



2 HabR 2 



dw n /dd 
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dw n /dR 

h 2 

2fx AB R 2 smH 



dw n /dip 



(64) 



As a matter of fact, the dominant contribution to V n (R c , R, if, t) is given by the first part of above 
equation (jHSJ). Magnitude of the second term (jSjj) can be generally considered to be small, due to 
presence of inverse mass factors M _1 and fi AB . Similarly as in the usual textbook treatments of 
the time independent Born-Oppenheimer separation 2^|, the correction (|64|) can be assigned to a 



difference between the Born-Oppenheimer and the adiabatic approximations. We shall neglect the 
term (J54"jl in our subsequent considerations, taking just 

d 

V n (R c ,R,i9,(p,t) = -ih — w n (R c ,R,$,tp,t) . (65) 
D. Electronic potential energy surfaces: The dc- field limit 

Let us investigate now a special case when the electromagnetic field varies only very slowly in 
time. Such a situation is referred as the dc (direct current) limit, and is encountered whenever an 
external electrostatic field <j)(f, t) is turned on and off sufficiently slowly, while the central frequency 
ujl of the laser pulse A(f, t) acquires a sufficiently small value (small with respect to the characteristic 
timescale of the electronic motions - i.e. practically even for microwaves). Theoretical analysis of 
this problem is facilitated by introducing a switching parameter rj (t) which modulates the slow field 
variations according to a formal prescription 

E{f,t) = E dc (r,r)(t)) , di](t)/dt^0 . (66) 

For the sake of clarity, we quote in this context a simple example of a low (microwave) frequency CW 
light E(f, t) = Eq cos(LVit — ki ■ r) with \ki\ = ujl/c and Eq ■ ki = . Here the switching parameter 
r)(t) can be defined just as uj^t and is then interpreted as an indicator of an instantaneous value of 
the electric field strength. 

Taking into account the property ()66|) of the field, it is clear that also the Hamiltonian interaction 
term (jlHj) depends on time only through rj (t) , and can be thus denoted by an extra symbol 



Wf(R c ,R,^<p, V (t)) = W el (R Cl R,&,<p,t) = - [M(&,(p)D%(R,0,<p,T«)\ ■ E dc (R c ,rj(t)) . (67) 
If so, the dynamical electronic states (|45J) take an explicit functional form predicted by the well 



known adiabatic theorem 



2l|. It holds 



$n(? > t; R c , R, if) = exp 



% - f £: d (R c ,R^^, V (r))dr 

a Jt a 



where the ?7-adiabatic electronic energies and wavefunctions are defined by an eigenvalue problem 



H^+WfiR^R^^^^l^ir^R^R^^^) = £^(R C , R,#,tp,v) , Rc, R, <P,v) ■ (69) 

Since the electronic Hamiltonian of equation ()69|) is real, the associated ?7-adiabatic electronic eigen- 
states $f (r N ; R c , R, if, rj) can be also considered as real quantities. Therefore, the matrix element 

(^\{d/dR)\$f) iN = ; (70) 



and, consequently, 

(<S> n \(d/dR)\$ n )? N = -^J^S^(R c ,R^,^r)(r))dr . (71) 

Analogical relations are valid also for the other matrix elements appearing on right hand sides of 
Eqs. (J59j) -([S0 |) . This shows that an appropriate solution of the problem (J59|) - (J6T3j) possesses the form 

w n (R c , R, 0, if, t) = l - f £f(R c , R, <p, 7/(7-)) dr ; (72) 

and the corresponding Born-Oppenheimer electronic potential energy surface reads as 

V n (R c ,R,ti,<p,t) = S^(R c ,R,#,<p,v(t)) ■ (73) 

Equation (f?3*j) displays of course an intuitively expected result: The calculated electronic potentials 
coincide with those obtained within the conventional Born-Oppenheimer approximation for the time 
independent static field. 



E. Electronic potential energy surfaces: The ac-field limit 

Here we investigate another special case when the electromagnetic field oscillates rapidly and 
(quasi)periodically in time. Mentioned situation is referred as the ac (alternating current) limit, and 
is encountered whenever the studied molecule is exposed to an adiabatically switched continuous 
wave UV-VIS-NIR laser, 

E(r,t) = Ejf,r ) (t))e + ^ t /2 + E*Jr, v (t)) e*>*/2 . (74) 

Here, symbol uo^ stands for the laser frequency, and the field amplitude E^ c (f,rj(t)) is allowed to 
depend very slowly (adiabatically) on a formally introduced switching parameter r\ (t) . To avoid 
confusion, we quote in this context a simple example of a CW-like Gaussian laser pulse E ac (f, rj(t)) = 
Eq e~ at2 e~ %kh ' r with real parameter a — > +0, \ki\ = ^l/c and Eq ■ ki = 0. Here the switching 
parameter r)(t) can be set to e _cr * 2 and is then interpreted as an indicator of an instantaneous value 
of an envelope of the considered light pulse. For the sake of completeness, we note also that the 
electrostatic potential 4>(r, t) has been chosen to be zero in the present example. 

Taking into account the property (|74|) of the field, it is clear that the corresponding Hamiltonian 
interaction term depends on time only through e ±luJLt and r)(t), and can be thus denoted as 



Wl(R c ,R,$,(p, v (t),t) = W el (R c ,R,$,<p,t) 



M{4,<p)D%{R,0,<p,i N ) ■ EJ?,r,{t)) /2 + c.c. (75) 



If so, the dynamical electronic wavefunctions (|45|) take an explicit functional form predicted by the 



adiabatic theorem for the Floquet states 



ljj. It holds 



$ n (r N ,t;R c ,R,fi,ip) = exp 



i 

h 



f S*{R c ,R,#,ip,r]{T))dT 

J to 



& n (r»;R Cj R^^, V (t),t) ; (76) 



where the ?7-adiabatic Floquet quasienergies and eigenfunctions are defined by a generalized eigen- 
value problem 

H el (i?) + Wl(R c ,R,$,<p,V,t) - ih{d/dt)] <S>Z(r N ;R c ,R,#,ip,r),t) = 

= £*(R C , R, tf, ip, V ) <(r^ R c , R, ^ V , t) . (77) 

In above equation (J77J), the time variable t is treated as an additional dynamical coordinate subjected 
to a boundary condition 

$£(r*i2 CJJ R, tf,<p, 77, t) = ^ F n (r N ;R cl R^,ip lV ,t + T) , T = 2ir/u L ; (78) 



in accordance with the spirit of the Floquet and (t,t') theories 24|. The Floquet wavefunctions can 



be expanded using the field free electronic basis set (|4*7jl into a sum 

m=+oo 

<S> F n (r N ; R c , R, 0, tp, V ,t) = £ E <w$» R, *, V, v) *U**, A) e +TOi ; (79) 

n' m=— oo 

which is reduced just to a single term ^(r^; i?) as soon as the field amplitude is turned off. Written 
mathematically, 

C^ nl (R c ,R,#,cp, m ) = 6 nn ,6 m0 for Ejf, Vo ) = . (80) 

Similarly, also 

£*(R c ,R,#,<p,r]o) = £n ■ (81) 

Strictly speaking, the above outlined formulation of the Floquet theory is physically adequate only 
in the weak field regime where the associated Floquet states resemble the properties of the bound 
states. For strong fields, where the field induced ionization phenomenon becomes important, the 
problem must be addressed within the framework of the non-Hermitian quantum mechanics, with 
different types of complex scaling transformations being employed to yield complex quasienergies 
(25 1 . In such a case an imaginary part of the quasienergy corresponds to an inverse lifetime of the 
associated metastable electronic Floquet state. Further details regarding the Floquet theory can be 
found in Refs. Q- 

Using an adiabatic ansatz (J76)) we find that the matrix element 

(* n \{d/dR)\* n )*< = -'-^J^SliR^R^^dr + (& n \{d/dR)\& n )iH . (82) 



Analogical relations are valid also for the other matrix elements appearing on right hand sides of 
equations ([5§ j) -([6Tl |) . Since the system of differential equations (JoT^ - lp^) is linear and homogeneous, 
an appropriate solution should possess the form 

w n {R c ,R,$,^,t) = \ [ E*(R c ,R,$,(p,r)(T)) dr + w n (R c ,R,$,<p,t) ; (83) 

where the quantity w n (R c , R, (p, t) is assumed to satisfy a system of the first order partial differential 
equations 

w n (R c , R, 0, <p, t) = - KI^rJO^ , dwJdR = -(<S>Z\{d/dR)\$l)? N (84) 

dw n /d# = -($*\(d/M)\K)t» , dw n /d<p = -($l\{d/dy)\<S>l) iN . (85) 

The question on solvability of Eqs. (jHH)-(jHSl) does not seem to be less difficult than in the case of 
Eqs. ()59 |) -()60 |) . For this reason, we prefer to carry out the phase transformation (|58J) using only the 
factor 

w n (Rc, R, <p, t) = l - f £*(R C , R, 0, <p, V (t)) dr . (86) 

Since the function (jHfij) does not represent an exact solution of the problem (j59jl - (j6()jl . the first 
derivatives of the hamiltonian term h^(t) (see Eq. (J53j) for n = n') are not completely eliminated. 
Instead, hW(t) is transformed into 

= - jj KlVaK)** • V^ c - £- (<S> F n \(d/dR)\$Z)? N A (87) 



VabR 2 n n dd fi AB R 2 sm 2 $ nlv ' Yn n/r dip ' 

Due to presence of the above first derivative term in the effective nuclear hamiltonian, the obtained 
electronic potential energy surface 

V n (R c ,R,#,ip,t) = £*(R c ,R,#,ip,ri(t)) (88) 

accounts only partially for the underlying Born-Oppenheimer nuclear dynamics. Nevertheless, as 
being explained in the next paragraph, the transformation (}58|) with the phase factor (}86|) still 
proves to be an important step which offers a lot of physical insight. 

The nature of the Floquet wavefunctions (J79*|) reveals that all the diagonal matrix elements con- 
tained in the formula (|57jl depend on time solely through the oscillating factors e ±imuJLt (m nonzero 
integer). Validity of this statement can be most directly demonstrated through an evaluation of the 
time averages 

rt+T 

J ($P(T)\(d/an)\$*(T))i* dr = o , n = (x,y,z, r,#, v ) ■ (89) 



by taking advantage of the Feynman-Hellman theory adapted for the Floquet states |18[. Due to 
presence of the just discussed oscillatory factors, the quantity (J88J) can be regarded as a physically 
well justified Born-Oppenheimer potential term in the ac-field limit. This holds true as long as the 
duration T = 2tt/ujl of one optical cycle remains much smaller than the characteristic timescales of 
the nuclear motions, so that the oscillating first derivative contribution (|KTj) is irrelevant. On the 
other hand, the correction ()87|) becomes increasingly important as uo^ decreases and approaches the 
dc- field limit. One can see it immediately also from the fact that the above ac-field potential ()88|) 
depends on time only through an envelope of the light pulse, while the correct dc-field formula (J73|) 
is defined in terms of an instantaneous field strength and includes thus the field oscillations e ±%UL *. 
The coming Section IV provides a more explicit comparison of the Born-Oppenheimer potentials in 
the dc-field and ac-field limits. 



IV. DIATOMIC MOLECULE IN AN ELECTROMAGNETIC FIELD: 
ELECTRONIC POTENTIAL ENERGY SURFACES IN THE WEAK FIELD REGIME 

A. Application of the perturbation theory in the dc-field limit 

Provided that the used field strength is sufficiently weak, the corresponding Born-Oppenheimer 
electronic eigenenergies and wavefunctions defined by equation 1)69)1 are only slightly different from 
their field free counterparts. If so, the framework of perturbation theory offers a straightforward 
method for resolving the mentioned Born-Oppenheimer electronic problem, and leads towards ex- 
plicit analytic formulas for the field induced corrections of the associated potential energy surfaces, 
£^ d (R c , R, ip, rf) — £n(R) ■ Clearly, the field strength is considered here to be the perturbation ex- 
pansion parameter, while the solutions of the field free eigenproblem (|%7|) are taken as an unperturbed 
reference. 

Before presenting the details of the perturbation approach, we would like to emphasize that this 
is not the only practical method for evaluation of the desired dc-field Born-Oppenheimer electronic 
potential energy surfaces. An alternative possibility is to use direct (ab initio) numerical solution 
of equation (|fi9)l. which is of course much more demanding from computational point of view, but 
remains appropriate even in the case of strong dc-fields. 

Perturbation expansion for the energy eigenvalue of equation (jo^j) possesses the form 

S^(R c ,R,^<f,r,) = £° n {R) + £^\R C , R, 0, <p, rj) + £^' 2 (R C , R, 0, <p, rf) + •■■ ; (90) 
where the first order correction 

£f'\R c ,R,$,ip,rj) = (<S> n \Wt?(R c ,R^,^v)\<)r» ; (91) 



i n 



and the second order term 



[$ n \W^R c ,R^^,V)\®°n)i 



£°(R) - £°,(R) 

Substitution of an expression (|dT|) yields consequently more explicit results 



(92) 



£f>\R c ,R^,<p, V ) 



M(#,tp)K\D*%(R,#,<p,?")\K)?» -E dc (R c , V ) ; 



(93) 



and 



£?' 2 (R c ,R,0,<p,ri) 



E 



(94) 



£°{R) - £« n ,{R) 

Relation ()93|1 determines of course an energy of an interaction between the molecular permanent 
dipole moment (associated with given electronic state n) and the total electric field (arising both 
due to the scalar and the vector potentials). We refer in this context to our discussion of the dipole 
term -D^g carried out at the end of subsection II. C. In passing we note that the just mentioned kind 
of dipole interaction constitutes a conceptual basis of the traditional rovibrational spectroscopy (we 
recall that the dc- field limit is appropriate in the microwave and the far infrared spectral domains). 
Physical contents of equation ()94j1 becomes apparent after rewriting it into an equivalent fashion 



£^ 2 (R c ,R,#,p,v) 
where the symbol 



M{ti,ip)El{R c ,r]) -a n {R)- M T ($, if) E dc (R c , rj) 



(95) 



a n (R) 



2 E 



£°,{R)-£°{R) 



(96) 



is recognized as a tensor of static polarizability associated with the n-th electronic state of the 
AB molecule. Since the quantity a n (R) is expressed with respect to the body fixed frame defined 
within subsection II. E, its off-diagonal matrix elements necessarily vanish provided that the relevant 
electronic state |$°) is of E type. (This is the case of the ground electronic states of the vast majority 
of diatomic molecules.) The situation is a bit more complicated in the case when the electronic states 
of types II, A, etc. enter into the play. We shall omit all the details for the sake of simplicity, and 
consider from now on just the E type electronic ground state |$°) . Symmetry considerations dictate 
that 



of (fl) = aWR) = a£{R) , a''(R) = aj(fl) , d?(R) = a* z {R) = af(R) 







(97) 



Moreover, it is clear that the x and y components of the dipole moment matrix element 
(<&g\D%g(R, i9, <p, r )|$°) ? iv appearing in the first order correction (J93j) are zero. In order to fully 



exploit these useful properties, we assume without any loss of generality that the only nonzero com- 
ponent of the dc electric field is pointing along the space fixed z-axis, i.e. 



E dc (R c , V ) = 0,0,E dc (R c , v ) 



Using equations (J2SJ), (|H7j) and (J9"5|) . we simplify the formulas (|9"3j) and into 



(9? 



£f>\R c ,R,#,tp, V ) 



E dc (R c ,rj) cos^($;| J D^ 2 (/?,^^,r 7V )|$») f 



(99) 



and 



1 



- El c (R c , rj) oMR) cos 2 $ + a^(-R) sin 2 ■& 



(100) 



Having derived these two important expressions, we conclude our discussion by writing down the 
final result for the associated translational/rotational/vibrational Hamiltonian (Jfi2|) assigned to the 
ground electronic state diatomic molecule *4£> interacting with a weak external dc electromagnetic 
field. It holds 



hf(t) 



h 2 d 2 



+ 



f 2 
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E dc (R c , n(t)) co S #(<S> g \D*% z (R, tp, r")\& g ) 



1 



- - E 2 c (R Cl r](t)) [al(R) cos 2 ■& + af(R) sin 2 ?? 



(101) 



Note that the dc-field enters here through its instantaneous time dependent strength E dc (R c ,rj(t)) . 
As being shown in the next subsection, this finding is in a sharp contrast to the situation when an 
ac-field is present. 



B. Application of the perturbation theory in the ac-field limit 

The framework of the perturbation theory can be exploited as well in the case of weak ac-fields. 
Perturbative treatment of the Born-Oppenheimer-Floquet electronic problem (J77J) leads towards 
explicit analytic formulas for the field induced corrections of the associated potential energy surfaces, 
S^(R C , R, i?, <f, v) ~~ £n(R) ■ Similarly as in the above dc-field case, the field strength is considered 
to be the perturbation expansion parameter, while the solutions of the field free eigenproblem (|47|) 
are taken as an unperturbed reference. In passing we note that, since the time variable t is within 
the Floquet formalism treated as an additional dynamical coordinate, we can still rely here on a 
(properly modified) time independent perturbation theory. Instead of supplying a more detailed 
description of such a Floquet type time independent perturbation approach, we refer to Q and H 
for well elaborated examples. 

Similarly as in above subsection IV. A, we would like to emphasize that the perturbation approach 
is not the only practical method for evaluation of the desired ac-field Born-Oppenheimer electronic 

•VI 



potential energy surfaces. An alternative possibility is to use direct (ab initio) numerical solution of 
the Floquet eigenproblem (|77jl . Such a task is of course much more demanding from the computa- 
tional point of view, but remains appropriate even in the case of strong ac-fields. 

Perturbation expansion for the quasienergy eigenvalue of equation (|77|) possesses the form 



F,2/ 



(102) 



where the first order correction £^ ,1 {R C , R, d, tp, 77) vanishes due to time averaging over the ac field 
oscillations, while the second order term 

2 



S*'\R c ,R,#,tp, V ) = £ 



N 



X 



1 
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Ei(R c ,v)/2 

1 



(103) 



£°(R) - S° n ,(R) + hu L £0( R ) - S° n ,(R) - Ku L j 
Physical contents of equation (jlOHj) becomes apparent after rewriting it into an equivalent fashion 



where the quantity 



M{#,tp)EL{Rc,ri)/2 -a n (R,u L )- M (1?, <p) -E ac (-R c , rj)/2 ; (104) 



a n (R,u L ) = - ^ 3 - - 2 -^-j- 



(105) 



is recognized as a tensor of dynamical (frequency dependent) polarizability associated with the n-th 
electronic state of the AB molecule, and an auxiliary symbol 



u nn ,{R) = \£° n (R)-£°,(R)}/h 



(106) 



Since the tensor a n (R,ui,) is expressed with respect to the body fixed frame, its off-diagonal matrix 
elements necessarily vanish provided that we restrict ourselves on the case of the ground electronic 
state |$°) = 1$°) which is assumed to be of E type. Symmetry considerations dictate then that 



a x g x (R,u L ) = af(R,u L ) = aj(R,u L ) , a?(R,u L ) = a\{Ruo L ) ; (107) 



and 



af{R,u L ) = a$*(R,u L ) = af{R,u L ) = . (108) 

In order to fully exploit these useful properties, we assume without any loss of generality that the 
only nonzero component of the considered ac electric field is pointing along the space fixed z-axis, 
i.e., 



EjR c ,r]) = 0,0, E ac {R c ,v) 



(109) 



Using Eqs. (J2B)), (fTTTTj) and (fTUHJ) . we simplify the formula (fTU^ into 

S*> 2 (R c ,R,$,ip,r]) = --E 2 ac (R c ,r]) [a\{R,uo L ) cos 2 ?? + o^,^) sin 2 ??] . (110) 

Note that we encounter here a multiplicative factor (1/4) contrary to the dc-field expression (jlOOj) 
where a factor (1/2) appears instead. Let us write down now the final result for the associated trans- 
lational/rotational/vibrational Hamiltonian (jfi2)l assigned to the ground electronic state diatomic 
molecule AB interacting with a weak external ac electromagnetic field. It holds 



^ {t) ~ 2M R < 2ii AB dR* + 2^ AB R 2 1 J 

1 



- -E^(R c ,r](t)) [al(R,u L ) cos 2 ?? + af(R,u L ) sin 2 tf 

Note that the ac-field enters here solely through its amplitude E ac (R c ,r](t)) , while the rapid field 
oscillations e ±luJLt do not appear. This behavior is in a sharp contrast to the case of dc-field analyzed 
within the previous subsection, see equation (jlOlj) . 

Our discussion concludes by pointing out one additional remark. The perturbational analysis 
performed within this subsection has been based upon the length gauge Hamiltonian defined by 
formulas (|41j ) - (j4~H|) . Another option might be to develop a perturbational expansion employing the 
momentum gauge Hamiltonian (fl"9|) . of course after having converted it properly into the body fixed 
frame following similar approach as in above subsections II.D-II.E. Interestingly, the results obtained 
within the length gauge and the momentum gauge perturbation theory possess somewhat different 
functional forms. Our calculations reveal that the momentum gauge expression for S^ ,2 (R C , R, t?, (p, t]) 
is still given by relation (|104J) . with the dynamical polarizability tensor being however redefined as 

The only difference between formulas (jl()5|) and ()112j) consists in the presence of an extra factor 
(u nn i(R) /ojl) 2 ■ Due to this extra factor, one would be tempted to argue that equations (jl()5|) and 
()112|) provide different results for the polarizability. Even a direct numerical calculation using a finite 
basis set of the field free molecular electronic states {|$°/)} shows that the predictions of the formulas 
()105|) and ()112|) are different. Especially pronounced discrepancies are found in the case when ul is 
far off resonant from any electronic transition |$°) — > |$°') . Mentioned observations might seem to 
indicate inconsistencies, since the passage from the momentum gauge to the length gauge is facilitated 
through an unitary transformation, and the quantum mechanical perturbation theory is known to 
be gauge invariant. A detailed theoretical analysis of the above sketched problem is postponed into 
another paper (^tJ. Here we just state without proof that both expressions f!105|) and ()112|) turn 
out to be identical, provided that one properly accounts for a complete set of all the electronic states, 



a n (R,uj L ) - t 2^ — 2 2 /m ~2 — • ( U2 ) 



including also the highly excited electronic continuum. Use of a truncated basis set in equation (jll2j) 
is not always justified and can provide an absolutely misleading outcome. This example should serve 
as a warning against an uncautious use of an incomplete electronic basis set (e.g. within the widely 
considered two level approximation) when studying the field induced atomic/molecular properties in 
the momentum gauge. 

V. CONCLUDING REMARKS 

In the present paper, we have studied general problem of the translational/rotational/vibra- 
tional/electronic dynamics of a diatomic molecule exposed to an interaction with an arbitrary external 
electromagnetic field. We have derived an appropriate body fixed frame Hamiltonian 1)4 lj). and in- 
troduced the concept of the time dependent Born-Oppenheimer approximation. An interesting open 
question on general existence of the time dependent Born-Oppenheimer electronic potential energy 
surfaces has been raised. Finally, we have derived an effective translational/rotational/vibrational 
Hamiltonian (jlOlj) resp. (jlllj) of a ground electronic state diatomic molecule in a weak dc/ac field. 
Our entire derivation is based upon the first quantum mechanical principles and well defined approx- 
imations. 

The theory developed in this paper is believed to be of importance for a variety of specific ap- 
plications, like e.g. alignment/orientation of molecules by lasers, trapping of ultracold molecules in 
optical lattices, molecular optics and interferometry, rovibrational spectroscopy of molecules in the 
presence of intense laser light, or harmonic generation. Moreover, the above outlined approach can 
be extended in a relatively straightforward manner to the most general case of a polyatomic molecule 
interacting with laser light. 
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